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Abstract 

We present a so-called adaptive Ising model (AIM) to provide a unifying 
explanation for sensitivity and perfect adaptation in bacterial chemotactic 
signalling, based on coupling among receptor dimers. In an AIM, an external 
field, representing ligand binding, is randomly applied to a fraction of spins, 
representing the states of the receptor dimers, and there is a delayed negative 
feedback from the spin value on the local field. This model is solved in an 
adiabatic approach. If the feedback is slow and weak enough, as indeed in 
chemotactic signalling, the system evolves through quasi-equilibrium states 
and the "magnetization" , representing the signal, always attenuates towards 
zero and is always sensitive to a subsequent stimulus. 
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As a prototypical system of cell signalling, bacterial chemotaxis is extensively studied in 
biology. A bacterium, such as Escherichia coli, swims under the control of several flagellar 
motors. When the motors rotate counterclockwisely, the bacterium runs forward; when the 
motors rotate clockwisely, it tumbles. The motors switch between these two modes of ro- 
tation with the probability ratio determined, through a signalling pathway, by the change 
of the concentration of the ambient chemical. Therefore the bacterium performs a biased 
random walk towards higher concentration of an attractant or lower concentration of a 
repellent. The signalling pathway (Fig. 1) is as follows [|IJ. The chemoeffector molecule lig- 
ands bind to transmembrane receptor dimers, which are coupled by two proteins CheW to 
two cytoplasmic histidine kinases CheA. CheA is autophosphorylated with the rate greatly 
enhanced by the receptor, hence attractant binding causes the receptor dimer to undergo 
a conformational change which leads to a decrease of autophosphorylation rate of CheA. 
CheA transfers phosphorylation group to two regulators CheB and CheY. Phospho-CheY 
directly modulates the motors. On the other hand, Phospho-CheB mediates demethylation 
of the receptor, while another regulator CheR promotes methylation. Attractant binding 
also makes the receptor a better substrate for CheR. Since methylation increases the au- 
tophosphorylation rate of CheA, Phospho-CheB and CheR provide a negative feedback. The 
basic structure of the information loop is: the ligand binding is the input, the activity change 
of CheA is the output, and there is a negative feedback from the output on the input. 

A crucial feature of chemotaxis is its sensitivity: as little as a single molecule can trigger 
a detectable motor response [[| . Another crucial feature is adaptation: after an initial sensi- 
tive response, the tumbling rate returns to the pre-stimulus level. A noteworthy fact is that 
there are about 2000 chemotactic receptor dimers clustering at a pole of the cell, furnish- 
ing a detecting "nose". Is the sensitivity enhanced by clustering 0? The possibility that 
ligand binding of one receptor can change the activities of more receptors was considered 
Q. A biological principle may be formulated: An attribute that exists most probably confers 
advantages over possible alternatives, especially if the latter have some apparent merit || . 
This principle and various experimental findings led to a cooperative model based on cou- 
pling among receptor dimers || . This model is equivalent to an Ising model in a bimodally 
distributed field, and provides an arbitrarily sensitive initial response, by choosing an ap- 
propriate value of a parameter comparing the coupling with the noise. In this theoretical 
framework, adaptation is achieved by a counteracting effect due to the negative feedback 
and mapping to an opposite induced field in the Ising model. The theoretical results are in 
good agreement with, say, a recent experiment ||. 

More investigations need to be made on adaptation and its relation with sensitivity. 
Especially it is important to explain why the adaptation is always perfect, i.e. the activity 
always returns to the pre-stimulus level precisely |7|]; recent experiments showed that this 
perfectness of adaptation is robust though other properties, such as the time needed to 
complete the adaptation, vary with conditions 0. 

In this letter, to improve the previous approach, we present a so-called adaptive Ising 
model (AIM), in which there is a negative feedback from the magnetization on the field. With 
large separation of time scales, there exists quasi-equilibrium, which is temporally local, i.e. 
on a short time scale. On a long time scale, however, the system evolves, with the existence 
of a dynamical attractor corresponding to a fixed pre-stimulus activity and is sensitive to 
the subsequent stimulus. Thus we obtain a natural explanation for why the adaptation is 
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always perfect, and show that not only the sensitive signal, but also the effective adaptation, 
is very likely a manifestation of the coupling among receptor dimers. Therefore sensitivity 
and adaptation are two sides of the same coin. Our model is above the level of molecular 
details, hence it is generic for signalling networks with similar conditions. Furthermore, the 
feedback from the information output on the input is a general way to preserve sensitivity. 

Under the assumption of high gain limit, the state of the receptor dimer, chracter- 
ized as Vi, can be either of two values, V° and V 1 , corresponding to the higher and lower 
rates of CheA autophosphorylation, respectively. As in the neural network ||, we assume 
McCulloch-Pitts behavior for Vi, i.e., in the absence of noise, Vi=ip(J2j TijVj + Hi — Ui), where 
ijj(x) = V 1 if x > while ip(x) = V° if x < 0, C/j is a threshold value, Ty describes coupling 
between neighbouring receptor dimers, Hi is the effect of ligand binding and methylation 
level change. With = Tji and Ta = 0, there exists a Hamiltonian which determines the 
equilibrium state. One may use the spin representation = 2(T^ — V°)/AV — 1, where 
AV = V 1 - V°. Defining J = J tj = 7^AF 2 /4 and B { = H { AV/2, assuming that the 
"magnetization" is zero at the paramagnetic phase for c = 0, and taking into account the 
time-dependence of Bi, we obtain the Hamiltonian 

= - E JaSiSj - E Bi{t)Si, (i) 

(ij) * 

where (ij) denotes nearest neighbouring pairs. = J > is a constant. An essential 
element of AIM is a negative feedback on Bi. 

— —— = -crbi{t-t r ), (2) 

where a > 0, t r is the retard time of feedback. This is not an arbitrary assumption put in 
by hand, but is a close representation of the experimental finding that, in the chemotactic 
signalling pathway, the state of CheA, through CheB and CheR, causes an opposite effect 
on its state later on |1J]. The initial condition is that Bi(t ) is bimodally distributed between 
B and 0, 

p[Bi(t )] = c5[B t {t = t ) -B} + (1- c)5[Bi(t )]. (3) 

Here Bi = B if the receptor dimer i is bound to ligand, otherwise Bi = 0. B > for 
attractant binding while B < for repellent binding, c is the the fraction of the receptor 
dimers with ligand bound, determined by the ambient concentration of the chemical. In 
other words, B^t) is superposed by two parts. One part is the externally applied field 
Bi(t )9(t — 1 ), where 6{x) = if x < while 9(x) — 1 if x > 0. Another part is an induced 
field, denoted as Mi(t), with dMi(t)/dt = —aSi(t — t r ). Note that we have assumed that 
the randomness is quenched, representing that ligand binding is so strong that debinding 
happens on a very long time scale f[U |. 



Generally, AIM defines a non-equilibrium model. However, the large separation of time 
scales holds in the current problem: ligand binding and conformation change occur within 
only millisecond, demethylation reactions take about 0.1 seconds, while time needed to com- 
plete the adaptation, which is associated with the slow modulation of methylation level, is 
on the scale of minutes [3J. This situation validates an adiabatic approximation. Roughly 
speaking, the demethylation reaction time plus the time for the transfer of phosphorylation 
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group from CheA to CheB corresponds to the retard time of feedback t r . Therefore, com- 
pared with the time needed to achieve temporally local equilibrium, the evolution of Bi(t) 
is very slow, i.e. the time scale of the adaptation process, characterized by t r , is very long, 
Hence we may solve our problem by using coarse-graining, replacing the above Hamiltonian 
with a temporally coarse grained one 

W(r) = - £ J tJ S t S, - £ Bi{r)S h (4) 

where r is the coarse grained and discretized time defined as r = int(t/t r ). Here the 
function int(x) is the greatest integer less than or equal to x. H(r) determines, through 
equilibrium statistical mechanics, the coarse grained instantaneous state characterized by 
the magnetization per spin, m(r), which is the value of each Si(r). Note that Si(r) = m(r) 
is the average of Si(t) over the time period from (r — 1)T to tT, equal to the ensemble 
average (Si)(r). On the long time scale, m(r) depends on r because of the feedback on Bi. 
m(r) is a measure of the signal. Usually the signal is characterized by the change of the 
number of the receptor dimers with state V°, the average of which is m/2. 

On the coarse grained time scale, the initial condition becomes that B^tq) bimodally 
distributes between and B, i.e. p[Bi(To)]=c6[Bi(r = r ) — B] + (1 — c)5[Bi(ro)], where r = 
int(t /t r ). The feedback equation becomes B^r) = B^t — 1) — ora{j — 1), or equivalently, 
M(r) = Mir — 1) — om(r — 1), which implies 

M(r) = M(r )-aJ2m(k). (5) 

TO 

On the coarse grained time scale, the induced field is the same for different spins, therefore 
the subscript i has been omitted. 

Under the adiabatic approximation, we apply mean field theory for each instant r to 
obtain: 

/ \ 2c 2 (!-c) , , cS 

mir) = I - - 1 (6) 

V; 1 + exp[-2/?(i/Jm(r) + M(r) + B)} 1 + exp[-2 ( 5(z/Jm(r) + M(r))] ' V; 

where M(r) is given by (ID, f3 — 1/ksT, v is the number of nearest neighbors. 

One may observe that m = is a fixed point of Eq. (|G|): if m(r — 1) = 0, then 
m(r) = m(r — 1) = 0. Moreover, if ojvJ is small enough, m(r) does not change the 
sign while its magnitude decreases towards [IT]. Therefore m = is an attractor of the 



evolution of the magnetization. 

In the original Ising model with c = 0, there are two phases, ferromagnetic and paramag- 
netic, depending on (3vJ . For AIM, however, as an interesting consequence of the feedback, 
m(c = 0) is always zero: suppose m(c = 0) is nonzero initially, the feedback automati- 
cally causes it to attenuate to zero. Therefore, we always have m(r < r ) = 0, and thus 
M(tq) = 0. Consequently 



m(r > r ) 



2r 



l+exp[-2/3(Wm(T)-0(T-To-l)(T m(k)+B)] 

£~1 K — Tf) 



_ i— ?M i (7) 

l+exp[-2P(uJm(T)-9(T-T -l)a'yi_ m(k))] 
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Thus when a "field" is applied, i.e. ligands are bound to the receptor dimers, randomly 
but with a certain occupancy c, there is an initial change of magnetization from to m(ro), 
depending on c. This initial response can be arbitrarily sensitive, as seen from dm(To)/dc 
with c — > 0, given by Eq. (10) of Ref. ||. However, due to the negative feedback of the 
output (magnetization) on the input (field) at each spin, the magnetization always attenuates 
towards zero. Practically, the adaptation is completed when the difference between m(r) 
and zero is below the detectable threshold of the motors. 

Note that in Ref. ||, it had to set that (3vJ < 1, i.e. the system should be in the 
paramagnetic phase when there is no ligand binding. With the feedback naturally integrated 
to the model in an ab initial way, this constraint becomes unnecessary, and thus the model 
becomes more robust. 

To obtain some analytical sense, consider high temperature limit (3 — > 0. In this case, 
m(r ) = c/3B/(l—(3uJ) ||. A simple calculation based on Eq. ([/]) reveals that m(r + Ar) = 
[1— (3a /(I— /5^J)] Ar m(r ). When (3a < \—(3vJ,m{j) attenuates towards zero exponentially. 
For generic values of the parameters, the solution can only be obtained numerically, as shown 
in Fig. 2. Note that the effective parameters are (3vJ, (3B, [3a, and c. Comparing plots 
for different values of parameters, one can observe that the speed of attenuation of m(r) 
increases with (3a and with /3vJ, while decreases with c. It increases with /3B, but when 
(3B is large enough, m(r) becomes independent of the exact value of (3B, as indicated by 
the results for (3B = 1, 10 with (3vJ = 0.5. On a log-log scale (not shown), the plots are 
generally convex, indicating that the attenuation is in general more rapid than exponential 
decay, due to the larger (3vJ. 

After the adaptation is completed, if there is a further change in the chemoeffector 
concentration, thus the occupancy changes from c to c + d at Tq, then m(t > Tq) is given by 
Eq. (0) with c updated with c + d. Because ?ti(tq — 1) = 0, M(tq — 1) is given by 

+ . r 2 <L-,!> (8) 



1 + exp[-2/3(Af(<^ - 1) + B)} 1 + exp[-2/3M(^ - 1)] 
Hence, 

m(r > T '\ = 2(c+c') 

V - 0J l+e^l-2^Jm(T)-aJ2l4,m(k)+M(rl,-l)+B)] 

o 

+ 2(l-c-c>) j (Q) 

l+exp[-2/3(y JmM-o-^L^ rn(k)+M(^-l)] V 1 

which is largely determined by d since the effect of c is counteracted by M(tq). m(r > Tq) 
attenuates towards zero, repeating the dynamics of Eq. (0). dm{r^)/dd with d — > 0, 
approximately equal to dm(ro)/dc with c — > 0, can be arbitrarily large if the latter can. 
Therefore our adaptation mechanism not only brings the signal to the pre-stimulus level, 
but also preserves the sensitivity, as required by chemotaxis. 

Therefore we have explained why perfect adaptation can always be achieved in chemo- 
taxis: a fixed pre-stimulus activity is a dynamical attractor. The variation of the values 
of the parameters, under a basic requirement that a is sufficiently small, only affect the 
time needed to achieve perfect adaptation. Thus our result is fully consistent with the 
experiments. 

Recent experimental analyses of the aspartate receptor revealed that attractant binding 
induces a displacement of one of four helices, each two of which constitute a subunit of 
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a receptor dimer Therefore Vj may be identified as the position of the mobile helix 

|i"3| . V° is the original position of the helix, corresponding to the higher rate of CheA 



autophosphorylation. V 1 is down towards the cytoplasm, corresponding to the lower rate of 
CheA autophosphorylation. Thus H is the force generated by ligand binding. 2B = HAV 
is the shift of energy difference between the two conformations due to free energy exchange 
with the bound ligand, or the work done by the generated force. One may find that 4 J / AV 
is the force due to the activity change of one nearest neighbour. 2Mj(t) / AV is the force due 
to feedback, and thus should be opposite to the force generated by ligand binding. 

In the high temperature limit, when At = —ln2/ln[l — /3a /(l — /3vJ)], m(r + At) = 
ra(To)/2. Assuming 1//3 ~ 4pN ■ nm, f3vJ « 0.5 ||, and that the time needed to complete 
adaptation be 1 minute, i.e. At rs 600, we may estimate that a rs 0.002piV • nm. Because 
the formula is for high temperature limit, the real value of a is smaller for the assumed 
values of the parameter values. Experimentally, by measuring (3, uJ, B and the adaptation 
time, a can be determined. On the other hand, a can also be determined through o = 
-[Mi(t) (t )]/ jj Si(t-t r )dt = [M{t q )-M{t)\/ E;ri ro m(r'). By comparing the results 
obtained in different ways, the model may be tested or refined. 

Eq. (0) implies that the feedback is assumed to be local. This is because we preserve the 
assumption that there exists a feedback loop for each receptor dimer although we consider 
coupling between the states of neighbouring dimers. However, one may make a straightfor- 
ward extension to include the neighbouring states in the feedback equation, without changing 
the qualitative physics. Furthermore, this makes no change in the temporally coarse grained 
feedback equation, or in the sense of (spatial) renormalization group. Therefore the large 
separation of time scales, which validates coarse graining, makes the essential mechanism 
not so much dependent on the microscopic details. This is also an aspect of robustness. 

Finally, it is interesting to note that in the case of bacterial chemotaxis, the clustering 
of receptors exists prior to a stimulus, while in many other cell signalling processes, the 
receptors diffuse on the membrane and the clustering appears as a response to the stimulus. 
One may make a generalization of our model to a sort of combination of Ising and lattice 
gas models to address such will be described in a forthcoming paper. 

To summarize, an adaptive Ising model is proposed to combine cooperativity and negative 
feedback. It is applied to the receptor network of bacterial chemotactic signalling and 
explains the perfect adaptation as a dynamical attractor. Both the signal magnitude and 
the sensitivity of response are adapted. The large separation of time scales leads up to 
the solution by using adiabatic approximation. The change of parameter values, under a 
basic requirement that the feedback effect is sufficiently weak, only changes the time needed 
to complete adaptation, without affecting its perfectness. Hence the robustness of perfect 
adaptation is explained. This work shows that coupling among receptor dimers gives rise 
to a unifying description of both sensitivity and effective adaptation. We anticipate further 
experimental and theoretical investigations. Cooperativity in cell signalling is likely a new 
playground of statistical mechanics. Combining cooperativity and feedback, and preserving 
sensitivity, the idea of AIM may be useful for a variety of problems. 

I thank Tom Duke and Dennis Bray for discussions. 
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FIGURES 



FIG. 1. A schematic illustration of the chemotactic signalling pathway of one receptor dimer. 
Around 2000 receptor dimers constitute a network. 

FIG. 2. Attenuation of m(r), the solution of Eq. (|7|), for different values of parameters, r is 
the coarse grained time, to is set to 1. To compare the attenuation speed for different values of 
parameters, we plot m(r)/m(ro). The parameters (fivJ, /3B, /3a,c) for each plot are given on the 
right upside. 
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